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We show in detail the general relationship between the Schrodinger equation approach to calcu- 
lating the MSW effect and the quantum field theoretical S-matrix approach. We show the precise 
form a generic neutrino propagator must have to allow a physically meaningful "oscillation proba- 
bility" to be decoupled from neutrino production fluxes and detection cross-sections, and explicitly 
list the conditions — not realized in cases of current experimental interest — in which the field theory 
approach would be useful. 
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I. INTRODUCTION 



Recent results from the Super-Kamiokande experiment ]lj have confirmed, with high statistics, the reality of the 
solar H and atmospheric Q neutrino anomalies. In the case of the atmospheric neutrino anomaly, the Super- 
Kamiokande data can be interpreted as providing evidence that the anomaly is due to neutrino flavor mixing [|J, 
and "long baseline" accelerator and reactor neutrino experiments may be able to confirm this conclusion J^j. Data 
] from Super-Kamiokande and the Sudbury Neutrino Observatory (SNO) eventually should allow tests to determine if 
■ the solar neutrino anomaly is also due to neutrino oscillations |fj. These developments are bringing to a critical test 
the possibility suggested by earlier observations of solar and atmospheric neutrinos — and by a persistent signal in an 
accelerator neutrino experiment, LSND |Q — that neutrino flavor mixing may provide one of the first experimental 
Q\ ' windows on physics beyond the Standard Model. 

In a generic neutrino oscillation experiment the event rate for the detection of vp from a flux of v a from a point 
f~| • source at distance L takes a form like the following: 



dT af} = / dE q ( L2d ^ a dEi ) {P^ ) , (1) 



where the direction of the neutrino momentum q points from the source to the detector. The first factor in the 
integrand represents the flux of neutrinos of energy i? q from a process involving a charged lepton of flavor a, and 
the third factor is the cross section for neutrino detection via a process involving a charged lepton of flavor j3. These 
factors are computed by the standard techniques of quantum field theory (QFT) , with the approximation of massless 
neutrinos. 

The middle factor — the so called "oscillation probability" — is typically computed with a quantum mechanical model 
of the oscillation process. In this model the neutrino state inhabits a Hilbert space whose dimension is equal to the 
number of neutrino flavors. The Hilbert space is spanned by a mass basis and a flavor basis. These bases are connected 
by a mixing matrix, taken to be the same as that which connects neutrino "flavor fields" to neutrino "mass eigenstate 
fields" in a field theory Lagrangian. The Hamiltonian is simply the particle energy, and the phenomenon of neutrino 
flavor oscillations arises because the Hamiltonian is not diagonal in the flavor basis. It was pointed out by Wolfenstein 
H that the parameters of flavor oscillations are altered in the presence of matter because of the effective mass induced 
by neutrino forward scattering off the background. The effective mass, which contributes to the Hamiltonian in this 
quantum mechanical model, can be computed by employing the famous formula relating the index of refraction to 
the forward scattering amplitude, where this amplitude is computed from QFT using standard interactions. Since 
the effective mass due to the background is diagonal in the flavor basis, in the limit of slowly varying background 
density the total Hamiltonian is diagonal in a new basis, the "instantaneous mass" basis. Mikheyev and Smirnov 
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subsequently noted that a level crossing of the "instantaneous neutrino mass eigenstates in matter" occurs in a 
background of monotonically varying density. The resulting "MSW effect" constitutes a new mechanism of flavor 
transformation that allows a parameter space of solution to, for example, the solar neutrino problem, that is very 
different from that provided by vacuum neutrino oscillations. 

The means just described of computing an experimental event rate for neutrino flavor transformations has the 
virtue of simplicity. However, being somewhat schizophrenic in its amalgamation of quantum field theoretical 
and quantum mechanical methods, it is not surprising that studies have appeared in which the neutrino produc- 
tion/oscillation/detection is examined as a single process in the context of QFT, with the neutrinos being virtual 
particles po[-^5[. These studies identify the conditions for which the amplitude for this overall process factorizes. 
However, by not showing the complete relationship between this amplitude and the neutrino production flux and 
detection cross section, these studies lack a firm justification (other than recognition from the usual quantum me- 
chanical picture) for calling a particular factor an "oscillation amplitude. "[|] In addition, previous works employing 
this "scattering approach" to flavor mixing have only considered vacuum flavor oscillations; it is our purpose here to 
consider the MSW effect in the context of QFT. Our approach is complementary to the work of Ref. Jl6| where this 
is derived from QFT within the context of relativistic Wigner functions. 

In this scattering approach to neutrino oscillations, the derivation of the usual Schrodinger equations for the MSW 
effect in the presence of a spatially varying background is a rather trivial consequence of the virtual neutrinos going 
on shell, since in that case the multiparticle nature of QFT becomes irrelevant. This property was noted in Ref. lO] 
in the context of vacuum oscillations. If this property is invalid, then the usual quantum mechanical treatment must 
be modified and the quantum field theoretical treatment becomes useful. Hence, perhaps the most important point 
of this paper is that the usual quantum mechanical treatment should be modified suitably if one of the following 
conditions holds: the on-shell neutrino momentum is nonrelativistic, the neutrino production or detection vertices 
are non-chiral, the wave packets of the external production/detection particles are sensitive to momentum variations 
of the order of inverse source-detector distance, or the neutrino effective mass splittings (determined by the effective 
potential including the background matter contributions) are large compared to the spread in the momenta of the 
external particles. 

We begin in Sec. II with a discussion of general neutrino oscillations. While this ground has been partially 
explored previously, the discussion will serve to clarify some of the physics of the oscillation process and identify the 
extent to which the neutrino propagator determines the probability of neutrino oscillations. In particular, we do not 
assume the vacuum propagator at the outset; we find the precise form a generic propagator must take to allow a 
physically meaningful "oscillation probability" to be decoupled from neutrino production fluxes and detection cross- 
sections, and pinpoint the component that gives rise to the oscillation amplitude. In Sec. Ill, we discuss the effective 
Lagrangian approximation assumed in our framework; as specific examples we discuss and neutrino backgrounds, 
including an outline of how to obtain a self-consistent neutrino background in, for example, the supernova environment. 
Having identified the Green's function (or propagator) as that which determines the portion of the neutrino oscillation 
probability that is independent of the production/detection mechanisms, in Sec. IV we study the Green's function of an 
effective theory of neutrinos in a static, uniform background, finding a rich pole structure. Great simplification occurs 
in the relativistic limit, and we recover the same oscillation amplitude obtained with the usual quantum mechanical 
model. In Sec. V, we study the Green's function in a nonuniform background potential. Under appropriate conditions 
we recover the usual Schrodinger-type equation for the oscillation amplitude. Sec. VI contains concluding remarks. 
An Appendix contains a quick and easy derivation in QFT of fully normalized neutrino oscillation event rates in the 
form of Eq. (pi), which is justified by the more complete treatment in Sec. II. 



II. GENERAL v OSCILLATIONS 



For the sake of completeness and to establish the setting of our calculation, we give a general overview of the 
neutrino oscillation calculation in field theory, going beyond previous works ]l0|-|l5| in a couple of ways. While 
earlier studies demonstrated the factorization of the amplitude under suitable conditions — enabling identification of 
the factor called the "oscillation amplitude" in the standard picture — they do not go all the way to a fully normalized 



1 A partial connection is made in Ref. [JLlJ, where it is shown, for example, how the 1/L 2 flux factor and on shell momentum 
space neutrino spinors arise from the vacuum propagator. However, the overall event rate they arrive at by "heuristic consider- 
ation," to use their words, contains a normalization constant. The relationship of this normalization constant to the coordinate 
space external particle spinors they employ is left unspecified. 
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expression like that of Eq. (Q), leaving one without an unambiguous, physically meaningful reason to name this 
factor an "oscillation amplitude." In addition, previous works exploring the conditions under which the neutrino 
propagator determines the oscillation probability, independent of the details of neutrino production and detection, 
have only considered the vacuum propagator. As a prelude to studying the MSW effect, we seek to elucidate the form 
a generic neutrino propagator must have to make it possible to disentangle the neutrino oscillation probability from 
the production/detection mechanisms. 

In this section we consider generalities without committing to a particular model (Lagrangian), establishing the 
connection between a physically clear definition of oscillation probability [given by Eq. (|J)] and the neutrino propa- 
gator. For a simple choice of external particle wave packets, we will give a fairly general expression for the neutrino 
oscillation event rate including a "calculated" normalization. For illustration, in the Appendix we give an explicit 
calculation leading to a fully normalized event rate like Eq. (|l|) in the context of a semirealistic model Lagrangian. 

In field theory and in physical situations, we distinguish a given flavor of neutrinos by their interactions with charged 
leptons. If we calculate the probability amplitude for the process involving flavor a at the source of the neutrinos and 
flavor (3 at the detector of the neutrinos, we can calculate the probability of neutrino oscillations of a — ► (3. The flavor 
of each interaction can be distinguished by measurable, on shell, external particles (i.e., the charged leptons). Hence, 
calculating neutrino oscillations is equivalent to calculating a scattering event where a neutrino propagator connects 
two flavor distinguishing vertices with external particles coming from them. 

In calculating scattering quantities such as cross sections, we usually calculate the plane wave scattering S-matrix 

S({k},{Pi}) ~ 1 = (2^) 4 <5 4 (^2(-l) d 'Pi +E(- 1 ) S!fc ') iM > ( 2 ) 

where M. is the usual invariant amplitude calculated with Feynman diagrams in momentum space and si = di = 1 
for incoming particles and for outgoing particles (the grouping of the momenta will be explained shortly) . This is 
considered to be a good approximation in the case of calculating usual collider event rates since there the events of 
interest occur within a single volume element before the final state particles are detected, and the corrections arising 
from localization of the interactions usually are not important to the detection rate. However, because neutrino 
oscillations involve quantum interference effects over macroscopic distances which separate the production point and 
the detection point of the neutrinos, we must take more care to account for the localization of the interaction points 
to calculate the leading order observable quantity. 

This localization of the vertices necessarily requires that the incident and final states be spatially localized wave 
packets instead of plane wave states. (As we shall see, analyses which appear to use only plane wave external states 
while restricting spatial integrations in an ad hoc manner |l2] , |l5| actually have complicated wave packets buried 
beneath the surface.) Hence, the probability amplitude^ is a superposition of Eq. (||), 

n l+Fn Is + Fs 

-4 = / JJ [dpj] ipDjipj,^) Yl [dh] ipsi{k,ki) [S({k m },{p m }) - 1] , (3) 

3 i 

where [dpj] = d 3 Pj/ [(27r) 3 ^/2£ , p J , {k m } are the external momenta of the vertex at the production region centered 
about xs, and {p m } are the external momenta of the detection region centered about yjj. Here the set of parameters 
p and k m characterize the peak of the wave packets' distribution of momenta.^ We have also fixed the number of 
external particles to be Is incoming and Fs outgoing particles at the source vertex, and 1 incoming and Fn outgoing 
external particles at the detector vertex. 

Now, let us see how this is related to the usual quantum mechanical treatment. The procedure for calculating 
P Ua ^un in Eq- (0) in a quantum mechanical model was described in Sec. I. In field theory, P Va ^v e defined by Eq. (|l|) 
can be calculated by comparing Eq. (Q) with the event rate derived from Eq. (0) . We associate As with the amplitude 

for {Is} — > z/J at the source and Ad with the amplitude for D+ vp—* {Fjj} at the detector (we specialize to 



Our conventions for the metric, gamma matrices, and normalizations are the same as Ref. @. The wave packet normalization 
is J d 3 p (2-7r)~ 3 |i/;(p)| 2 = 1. A plane wave packet i/j(p,p') = ^=-5 3 (p — p'), where V is a volume factor, follows this 
normalization convention provided [<5 3 (p — p')] 2 i s interpreted as j^jjS 3 (p — p'). 

3 The Oth component of these "parameters" is taken to be on mass shell since we will choose the wave packets such that they 
behave like plane waves for large values of the spatial components of these "parameters." 
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one detector particle D), and choose plane wave packets for the source's final state (anti)neutrinos and the detector's 
initial state (anti)neutrinos: 



. Is+Fs 



As = ^=jy J IT ^(hM [S S ({k m },q) - 1] , 

+ F D 

n W ^Dj(Pj,P } ) [S D ({p m }, q) - 1] , (4) 



where g = (Eq, L-Eq) is the neutrino momentum, and L = (yo — xs)/|yo — xg| points from the source to the detector. 
In „4s and ^ we have implicitly assumed that the neutrinos are massless, because massive flavor eigenstates cannot 
be asymptotic states. Standard kinematics then yields the relationship 

\A\ 2 f dEqEq \A S \ 2 V \A D \*V 

"M-) (-) m ' v°J 



T J {2-KfLH vD T s i-J^i-J T D ' 

where V is the usual total volume factor associated with the phase space and normalization of plane wave packets; T, 
T$, and Tjj are the usual time factors associated with stationary wave packets; v v d is the M0ller speed (associated 
with the flux) between the detector particle and the neutrinos; and L = |y£> — x^|. 

As just indicated, we make the simplifying assumption of stationary wave packets. The main simplifying utility of 
this energy conservation approximation is to get rid of the neutrino momentum integral.^ We encode our assumption 
of stationarity by defining spatially smeared functions 

. Is+Fs 



j 

r 1+F D 

ffD(y, e^i^ 'a-" = J 1] [dpA^DifapJe'ZS-V^-ViMDiipihq). (6) 



We note that A4s and A4r> have the form (assuming V-A lepton currents) 

M S ({ki},q) =u-(q)P R M 1 ({k i }), M D ({ Pi },q) = M 2 {{pi})P L u-(q) {v osc), 

Ms({h},q) =M 2 ({ki})P L v+(q), M D {{p l },q)=v+{q)P R M 1 {{p l }) (P osc), (7) 

where M\ and M 2 are respectively column and row vectors in spinor space, and Pl and Pr are the left- and right- 
handed chiral projection operators. Following the conventions of fl7}| , we represent the spinors u and v as 

where er M = (l,tr), = (1, — <x), and er is the three-vector of Pauli matrices. Since these are spinors for massless 
particles, the spin index s is associated with the spin component along the momentum axis; specifically, we have 
£- = T)+, = -TT, with a- ■ q £±(q) = if^q). Putting Eqs. (|,§) into Eq. (§) yields 

(2tt)5 fc(-l)"^ + EC" 1 )*^) P '-j^y = 16 ^ 2 ^ Iff5(-q, {fcj)f 2 |fe(q, {^f 2 , (9) 
where we have defined 

~g D (q, {£.}) = J dhg D {y,{p^,q)e<^i { - l)dl ^) y , (10) 
and similarly for gs. Stationarity constrains |q| = i? q = — ( — 1 ) Si -Ekj • 



We refer the reader to Ref. (l4| and references therein for related discussions regarding coherence. 
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Let us turn our attention to A. Given that we have a neutrino propagator G in our amplitude, and assuming V-A 
lepton currents, we can write 

S({k i },{p j })-1 = I d 4 ye^ { ' 1)d ' Pl V J d^xe^-^'^z J J^Le*"'^ 'M a P L G(s) P R M u (11) 

where M\ and M.% are the same as in Eqs. (^), and s is the off-shell propagator momentum. The upper (lower) 
sign of =p in the exponential is for neutrino (antineutrino) oscillations; this arises from choosing x (y) to always 
correspond to the source (detector). That is, for neutrino oscillations of flavor a to flavor (3, the Green's function is 
iG 0a (y,x) = (T{^(y)D a {x)}) = i J e" ls <y^G^ a {s) (with T{} and () denoting a time-ordered product and 

vacuum expectation value respectively), while for antineutrino oscillations a — > /3, the labeling is iG a ^(x,y). 
Insert the identity 

g M r '{1^1 A = 1 (12) 

on both sides of the Green's function in equation Eq. (pT|), where as before q = {Eq, L_B q ), and we define £ = 
(Eq, L^(Eq) 2 — m 2 ) in which the parameter m? < E 2 , though its precise value is unimportant in this context. We 
note that since q is null, q ■ 7 = J2 S uS (q)u s (q) ( or 1 ' 7 = J2 S vS (l)v s (<l) , if one wishes to consider antineutrino 
oscillations). From the explicit form of u and v it is easy to see that 

PL (q ■ 7)(S • 7) = PL u-{q)u~{q) (£ • 7), P L {£, ■ 7 )(q • 7) = (f • 7) P R u+ (q)u+ (q), 

K ■ 7) (9 • 7) ^fl = K • 7) PR, (q ■ 7)(€ ' 7) ^fl = « + (9)S+(g) Fl (£ • 7)- (13) 

[The same relations hold for u ± ((j') replaced by v^(q).] Soon we will show the form that the Green's function must 
have, after localization by the source and detector, in order that the term with u~~ (or v + ) on both sides of G be 
the only one to contribute. If more than one spin contributes, we will not recover the usual quantum mechanical 
treatment without spins taken into account. In that case, working with the full scattering picture of Eq. (|^) is useful. 
Keeping only the term with the relevant spinor on either side of G, one can show that Eq. (|^) becomes 

A=-j d*yg D fr,{2.U)jZS-V d %-v J fa g s (x,{k.},q) c'E.C- 1 )"*.-* i J *L e V<v-*) p G (a) P, (14) 

where P = ^°u-(q)/(2E tl ) = ^°v+(q)/(2E cl ), P = P^°, and g s and g D are given by Eq. (|). 

In passing, we would like to remark that in Eq. (pT[), we can always write (for neutrino oscillations, for example) 

g S (xAk},q)e^ ( - irik <- x =iu-(q)M 1 ({k i })f S (x,{k i }) (15) 

(and similarly for gjj) where f$ is a scalar function of x. Then we obtain the form (again assuming the single-spin 
contribution is justified after spatial integration) 

A = J ^ar^tf/s^, {&})/!?(»,{£(}) / ^^^'^({k^s^p.}), (16) 

which is a common starting point of analysis in the literature as in Refs. JT2| , |l5[ . However, with arbitrarily chosen 
smearing functions /, it is difficult to assess what actual scattering question the amplitude is an answer to, because 
the smearing functions are not the wave functions of the in-out particles, but are the wave functions smeared over 
the matrix elements. As a consequence, the normalization is usually ignored in this approach. We will return to the 
normalization later in this section. 

Returning to Eq. (|l4|), integration over a; , y°, and s° gives an overall energy-conserving delta function and sets 
s° = Eq (for antineutrino oscillations, s° = —Eq). We also note that the chiral structure of P GP (as well as the 
original matrix element) picks out only the Glr block of the neutrino propagator, where Glr is the nonzero 2x2 
submatrix left by P^GPr. In addition, the localization of gs and go around X5 and yn respectively "clamps down" 
on the coordinate space Green's function. In particular, if the characteristic widths L$ and Lp of gs and go are 
much smaller than the source-detector distance L = \yjj — xg|, we note that if the "oscillation probability" is to be 
disentangled from the details of neutrino production and detection, the relevant portion of the Green's function for 
oscillations a — > (3 must take the form 
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Gl a R (s° = E v y, x ) = J-0- e -'(y-x) G ^ (s o = ^ g) 



(2tt) 3 

/ „\ gi-EqL-(y-x) 

-£7 q (l-o-.L) — ?HP a (E v y Dt xs) (v osc), 



g^( S ° = -^ q ,x,y)~+£; q (l- --L) e ^ qL(Y X) y0|X9 ) (pose), (17) 

V / 47r|yz3 - x s | 

where L = {yjj —xs)/\yD — xg| points from the source towards the detector, and the quantities 7? and H have only 
flavor indices. The factor Eq ^1 — er • L^j arises from the kinetic term in the Lagrangian, and takes this form due to 

the relativistic limit. Another key ingredient is the factor e l " Eq ' y_x ' , which is the leading phase factor in the relativistic 
limit coming from e ls '( y-x ' evaluated at the poles of Glr- In addition, l/|y — x| comes from the asymptotic expansion 



of the left hand side of Eq. (17) in the limit that |y — x| — > oo, and it can be considered to be the monopole term in 
a multipole expansion. We will discuss the validity of the factorization and the asymptotic expansion further below 
in momentum space. 

Before we talk about momentum space, let us give an example of Eq. (|l^) by considering the vacuum propagator. 
In that case, it is straightforward to show that (anticipating the relativistic limit) 

G LR (s°, x, y) = (s° + ict-V) [M- l G RR (s°, x, y)] , (18) 

r ,,_i^ , o vi«/3 e l l s °H x -yl ^ / m||x-y|\ 

[ ^rr^^y)] = - 4?r|x _ y| u °s u h ex p y~ i 2 |s°i J ' (19) 

where Grr is the nonzero 2x2 submatrix left by PrGPr, M is the mass matrix appearing in the Lagrangian, and the 
rrij are the mass eigenvalues. In Eq. ( |l9| ) we have made the flavor indices explicit; the relationship between the flavor 
fields and mass eigenstate fields is v a — ^ Uaiipi, where the U a i are elements of a unitary matrix. For |s°|L 1, 

Gf fl ( S °,x,y) ~ -^^J [s° - \s> ■ f(x,y)] £ E/^-exp (^ "'j^' j . (20) 

where f(x, y) = (x — y)/|x — y|. To apply Eq. ( |2p| ) to neutrino oscillations one takes s° — > Eq and x, y — > y,x. For 
antineutrino oscillations, s° — > — _E q and x, y — > x, y. After making these substitutions we will be integrating Eq. 
( pp| ) over localization functions of characteristic widths L$ and Ljj centered on x = xs and y = yo. This means that 
for Ls, Ljj <C L, we may replace x and y by xs and yo everywhere except the phase factors, in which we consider 
the first order variation, 

|x - y| ~ |x s - y D \ + L ■ [(y - y D ) - (x - x s )] 

= L-(y-x). (21) 

We see that a necessary mathematical condition (in addition to the relativistic assumption and EqL ^> 1) for 
Eqs. ( |l8fpi| ) to reduce to the form of Eqs. (jL71) is m| Ls.d/ i^Eq) <C 1 [which also implies the more familiar 
(rrij — mf) Lg i £)/(2£' q ) <C 1]. The physical basis of these conditions can also be inferred from Eqs. (JT^jSTJ). EqL ^> 1 
allows the propagating neutrino to become an on-shell relativistic particle, and also allows appreciable oscillation 
phase to build up over the source-detector distance, rrij Lg t ]j/(2Eq) <C 1 requires that no appreciable oscillation 
phase build up on length scales comparable to the width of the external particle wave packets. The necessity of 
these conditions for disentanglement of the flavor oscillations from the details of neutrino production and detection is 
evident. 



The origin of these conditions can also be understood in momentum space. First, rewrite Eq. (14) as 



xhn(s, {£.}, q)h s (-s, {fcj, q)P G ^,s° = T , P 



(22) 



where the functions hs and ho can be approximated to have no xg or yu dependence. This can easily be seen to be 
exactly true for the ideal case of isotropic smearing functions, e.g. g_o(y, {p }, q) — hjj(\y — yu |); Eq. (PHI) then yields 



G 



9i5 = e lu ' yD ti£>(\u\), where u = s l) di P;- Because the propagator will in general have poles corresponding to 

the mass of the physical neutrino states, the dominant contribution to the integral in the asymptotic limit LEq — > oo 
will be from a term that contains the integrand of Eq. ( p2| ) as a factor evaluated at the poles and stationary phase 
points (critical points). For the vacuum, the constant potential, and the adiabatically spatially varying background 
potential cases, one can asymptotically expand the integral Eq. ( [22] ) in the limit LEq — > oo (similarly as in Ref. [[TT||) 
to find that to leading approximation, the term /i_d(s, {p^}, q)hs{—s, {hi}, q) can be moved outside of the integral with 

the replacement s — > s* where s» corresponds to one of the critical points. By factoring out h, we have implicitly 
assumed that hu(s,{p.},q)hg(~s,{k i },q) is not sensitive to the splittings in the critical points (otherwise different 

pole momenta s of Glr will cause /i_d(s, {j^}, q)hs(—s, {fc;}, l) to have different values, preventing factorization). This 
means that the wave packets must be flat in momentum space at least within the range of pole momentum splitting. 
Also, if this is not the case, one of the poles will not contribute (because the amplitude of the wave packet has fallen 
off with respect to the amplitude at the other pole), and no neutrino oscillations will occur (or more accurately, the 
neutrino oscillations will be greatly suppressed relative to the background). We will refer to this flatness of the wave 
packet as insensitivity to s* splitting. 

Also note that because of the presence of the exponential in Eq. (18), this leading term in the asymptotic expansion 
will not be a good approximation unless the inverse "momentum scale height" (i.e. logarithmic derivative) of ho hs 
near the poles is much less than L. Hence, factoring out the wave packet dependence which is crucial for the validity of 
the usual quantum mechanical treatment requires the wave packet factor ho hs to be insensitive under 1/L momenta 
variations as well as the splitting variations. 

While localization is clearly necessary for the observation of oscillations, the source and detector localization scales 
L5, Ld implied by Eq. (^]) cannot be smaller than the Compton wavelength of the lightest external particles. In the 
case that all the external particles connected to a given vertex are nonrelativistic, this gives rise to a constraint on 
the masses of these external particles. To see this, consider the ideal case mentioned above in which gu(y, {p^}, q) ~ 

hody — Yd\)- Then go = e tu ' yD hn(\u\), where /i£)(|u|) is damped for |u| = |s» — J2i(~ l) di P ; l larger than 1/Lp. 
Hence the splitting insensitivity condition can be written as 



X)(-irk,+s. 

£(-1)^-3* 



< — < M LS 
Ls 



< — < M LD . 
Ld 



(23) 



where we have denoted the lightest external particle masses to be Mls and Mld for source and detector, respectively. 
The critical momentum will be s* ss L^/(y^,(— l)">Ek, ) 2 — rh 2 where rhj is the effective pole mass of the particle. 



Hence, if rhj < \ J2i(~ l ) s ' E k 



V^(-i) si ^) 2 - 

and the external particles are nonrelativistic, then Eq. 



t) can be satisfied only if 



about equal mass of external particles enter and leave the source/detector vertices. If any of the the external particles 
connected to a given vertex are sufficiently relativistic, this severe constraint does not arise. 

We now show that Green's function must take the form found in Eqs. (|l7} ) after being spatially "clamped" by the 
source and detector if the terms projected by the spinors u~ (or v + ) on both sides of G are to be the only contributions 
to the amplitude. In Sec. IV, where we study the neutrino propagator in a uniform, static medium, we will find 
it convenient to identify L of Eq. (|l^) with the positive third spatial direction. In that case it is straightforward 
to show, using Eq. (|l^), that the terms with spinors of different spins on either side of G pick out the off-diagonal 
spinor space elements of Glr, while the terms with the same spins on both sides pick out the diagonal spinor space 
elements. The matrix (l — cr 3 ) from Eqs. (17) confirms that only G\ 2 R is nonzero, and therefore only the term with 
vT (q) [v + (q)] on both sides of G survives for the neutrino (antineutrino) oscillations.^ 

Recalling that q = Li? q , and upon inserting Eqs. ( p7| ) into Eq. ( fl4| ) and Eq. (14) into Eq. (||), we finally arrive 
at the neutrino oscillation probability 



— I TlP a 



H^iE^yo^s) 



P^ e = |jj«0(E q ,y D ,xs)| 



[y osc), 
(y osc), 



(24) 



3 Note that with L set to the third spatial direction, both u (q) and v + (q) have 4-spinor components (0, 
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where H and H are defined by Eqs. (0), and we have assumed that the detector particle D is nonrelativistic 
such that v v d = 1- With the cancellation of the source and detector wave packets, one can see why employing 
a separate quantum mechanical model to compute the oscillation probability is possible. [Note that the standard 
vacuum oscillation probability is recovered here, as is clear from Eqs. ( |l7j - |2"c| )1 . 

We emphasize that the Green's function here is the full propagator in any given theory and we have made no 
severe assumptions about the nature of the production and detection effective vertex.^ Hence, as expected, the 
production/detection independent field theoretical effects on neutrino oscillations come from the coordinate space 
Green's function. What perhaps is less expected is the fact that unless the wave packets of the external particles 
satisfy specific properties, the transition probability will not just depend on the propagator, but the entire coherent 
scattering process which neutrino oscillation really is. Such tangled wave packet dependence is discussed for example 
in Ref. @. 

Before concluding this section, we would like to note that we can easily work out the neutrino oscillation detection 
rate including the normalization if we assume a particular class of wave packets. Let us define a box wave packet as 
a configuration such that the superposition integral gives, for each outgoing particle, for example,^ 



d 3 k 



--4>s(k,k)e 



iE k x° 



d 3 k 



(25) 



where N is a constant independent of x and B(z) is a function which vanishes if z is outside of a box centered about 
the origin with each dimension of length Lg and is 1 everywhere else. The approximation in Eq. ( p5| ) is valid for 

x° <C (£k£s)/(27r|k|). The normalization condition J -^^^{k)] 2 = 1 fixes N and implies 



Mk) = 4= 
v V 



1-0 



1 



I?E\ 



i(k-k)-x s£>s ( k _k) ) 



(26) 



where 



D s (v) = 8 



su\(v x L s /2) sm(v y L s /2) sin(v z L s /2) 

V X VyV Z 



(27) 



Hence, since the box scale must be larger than the Compton wavelength scale, these external particles will generally 
have "plane wave in a box" type of normalization [up to a (27r) 3 5 3 (k — k) type of localization factor Ds]- This wave 
packet can be used to calculate the event rate using Eq. (O) in a standard way. Since ipg will have a width 2tt/Ls, 
smearing of any function that is proportional to momenta whose magnitude at the peak of the distribution is of the 
order 2ir/Ls, (or less) will deviate significantly from the Dirac 6 smearing of that function. Fortunately, because Mi 
and M2 do not depend on such small momenta, we can write the amplitude in the form of Eq. (|l6|) with 



Is+F s 

f s (x,{ki})= n 



2E k .V s 



-S(x-xs) 



(28) 



(and similarly for fu) which is what one would use to calculate scattering of particles confined to a box interacting 
with particles that can propagate outside of the box. Explicitly, the transition rate per source and detector particle 
is given by 



dT = (2t:) 8 \^2(-irE^ 
d 3 s 



1 




1 



2E k V s 



F D +1 

n 



1 f| fMs ff d 3 P a V D 

11 (0^3 11 



2E P V D 



(2tt) 3 



3 ±is-(yr)-x s 



(2tt) 3 
2 



(2rrY 



(29) 



The most significant assumptions leading to our final result were the V-A type of lepton currents; the stationary approximation 
in Eq. (^|); relativistic neutrinos, i.e. (_E q — |s,|)/_B q <C 1, where |s»| is the magnitude of a pole in the momentum space 
propagator; sufficiently localized and separated source and detector, i.e. (_E q — IsjDLs^ -C 1; and E^L 3> 1. 
7 We will write the wave packet centered about xs, since the one centered about yzj is analogous. 
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where D$ is defined by Eq. (P7|), and Djj is similarly defined. In this case, one can also use the usual heuristic box 
quantization formalism to calculate the event rates including the normalization. For pedagogical purposes we carry 
out this simple exercise explicitly using a fermion field toy model in the Appendix. 

To summarize this section, we have shown to what extent the neutrino "oscillation probability" is determined by 
the production/detection wave packet-independent propagator of the field theory. If wave packets for the production 
and detection events described in Eq. ([!]) satisfy suitable localization properties and the effective mass splitting of 
the neutrinos is not large compared to the momentum width of the wave packets, the neutrino propagator determines 
the probability of transition as defined by Eq. ([!]). This factoring of the wave packets out of the transition amplitude 
is crucial to recover the usual quantum mechanical picture of neutrino oscillations. Furthermore, we see how the 
multiparticle nature of the field theory becomes irrelevant as the poles of the propagator are the only states to 
contribute in this limit. For this factorization to be possible, the source-detector separation L must be large enough 
such that the wave packets do not vary over l/L momentum perturbations about the pole momentum, and the pole 
momenta splitting must be small enough such that the wave packet amplitudes take on approximately the same value 
for the various pole momenta [as discussed between Eq. (|2^) and Eq. (plf)]. Furthermore, since the usual quantum 
mechanical treatment neglects the spin of the neutrinos, only one spin projection of the Green function must contribute 
to the amplitude to recover the usual treatment. We have seen in this section that the relativistic limit of the on 
shell neutrinos and the chiral nature of the interactions ensure this. Now that we see that wave packet dependence 
can be factored out (as is implicit in the usual simple quantum mechanical treatment), we shall concentrate on the 
wave packet-independent field theoretic calculation of the MSW effect, which is encoded in the propagator within a 
background medium. 



III. EFFECTIVE LAGRANGIAN 



In this section, we briefly explain the effective potential employed in our calculation of the MSW effect. Focusing 
on the physics well below the electroweak scale, we write the usual electroweak effective Hamiltonian density as (see 
for example |l8|]) 

-Hl=% (.#4, + AJn,) , (30) 

where is the charged current and J 1 ^ is the neutral current. Take for example the contribution to the neutrino- 
electron interaction of the form 

Hi = —^-D el »P L v e e ltl P L e, (31) 

which will be dominant for the MSW effect. We can distinguish two different types of scattering: forward scattering, 
for which the background particles do not change their momenta; and non-forward scattering. In calculating our 
transition rate, we will not account for non-forward scattering contributions because these can be considered to be 
separate production events. With this restriction, in expanding the S-matrix perturbatively, the main background 



contribution will come from the expectation values of Eq. (31) taken with respect to the electron background states 



The scattering amplitude will then receive contributions proportional to powers of 

(VP = 2V2G F (n\e ltl P L e\n) (32) 

where n labels some many-body background electron state (not necessarily translationally invariant). Note that the 
right hand side is proportional to the left handed electron current of state \n). In an experimental setting, we are 
really interested in ensemble averages of the probabilities (not the averages of the S matrix) . However, for macroscopic 
numbers of electrons, we expect the main contribution to come from a set of degenerate states having the same spatial 
localization as the macroscopic distribution function. This approximation will break down if the density matrix is not 
sharply peaked about one set of states giving degenerate contributions to the scattering amplitude. We will assume 
that such a peaked distribution exists, and we will merely assign macroscopic currents to the expectation value of 
currents that will arise in the scattering amplitude calculation. This means that we will replace the interaction 
Hamiltonian density of Eq. ( |3l| ) with the effective density 

nf = ^-v el »p LVe j; (33) 
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where is the macroscopic left-handed electron current. For example, for an unpolarized e ± background one would 
employ — based on consideration of the sum over spin states of single particle expectation values of e^P^e, for 
example — the following expression: 

V = kf$F\fr(P)-f*W]%> (34) 

where the / e ±(p) are the usual distribution functions, including a factor of two for spin degeneracy. As usual, this 
procedure neglects higher order correlations. Note that for electrons in thermal equilibrium, our prescription, e.g. 

(V°) =V2G F (n e - -n e+ )^ (35) 

gives the same mass shift as the real time thermal field formalism employed by |L9| . 

As another example, the effective potential due to background neutrinos is of interest in the envelope of a super- 
nova/nascent neutron star, where the neutrino flavor composition can affect, for example, the explosion mechanism 
pO| or the outcome of possible heavy element nucleosynthesis [21 22| . In addition to the e ± background, we must 
consider neutrino- neutrino forward scattering arising from another term in Eq. (^0|), 

Hi = — = ^ Vil^PLViVjluPLVj, (36) 

where the indices i,j label the mass eigenstate fields. We work in the mass basis because the external neutrino 
background consists of on shell states, a point whose consequences were emphasized in Ref. |^| (see also Ref. p2[ , 
and references in these) . As seen previously, in the perturbative expansion of the S-matrix we will have occasion to 
take a background expectation value of this interaction (this time with respect to a many-body background neutrino 
state). Two of the neutrino fields will be paired with fields in the "production" and "detection" interactions, leaving 
two other fields whose background expectation value is taken: 

< W '> = 7§ E 2D,YPlM^Pl^) + ^§ E ^PLiWjh^PLVj. (37) 



While the correspondence of the expectation value in the first term of Eq. (|37]) with a macroscopic current is apparent, 
the meaning of the second term is less clear. 

We seek guidance by considering the expectation values with respect to the single particle neutrino states \qsvk) 
of momentum q, spin s, and mass m^. The expectation value in the first term of Eq. ( |37|) is 

6 

(qsi/fcl^Pi^-lqsi/fc) = - f u{qsv k )^P L u{qsv k ). (38) 

We consider a relativistic neutrino b ackg round, so that to leading order there are only negative helicity states; then 
the momentum space spinors in Eq. ( p8[ ) are approximately 



where (6, (j>) denote the polar and azimuthal angles that define q. The first term in Eq. ( p7| ) becomes 

^^2^P L u i (D j ^P L , j )=V2G F ^x\^^Xi, (40) 

i,j i,j q 

where Xi denotes the upper two components of PiMi and = (1, — er). Turning to the second term in Eq. (|37j), one 
finds 



qsiy k \viVj\qsis k ) = - toJ\3{oj^ \ "(l s v k)u{qs v k ). (41) 



Si k Sj k 
(2^B q ) 

Employing Eq. (|39|), the second term in Eq. ( |37| ) becomes 
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Noting the similarity between Eqs. (|7J) and (fl2]), in the relativistic limit we replace the Hamiltonian density of Eq. 
( |36| ) with the effective density 

Hf = V2G F nJ^P L n + V2G F ViJ^P L Vi, (43) 

id i 

where 

J ^ = J^ifM-M(p)]^. (44) 

The flavor fields v a are related to the mass eigenstatc fields vi by v a — U a iUi, where the U a i are elements of a 
unitary matrix. In the limit of vanishing mixing angles (£7 a i = S a i) and a thermal background, the effective interaction 
of Eq. pg]) gives the same mass shifts as obtained in Ref. [^9|. For nontrivial mixing, however, in terms of the flavor 
fields Eq. (j43|) becomes 

Hf = V2G F ^J^l^P L y a + V2G F (EW£ U l(j ) ^P L v p . (45) 

ad a.,0 

Considered as a matrix in flavor space, the quantity in parentheses in the second term of Eq. ( [45| ) contains off-diagonal 
elements. It is clear from our derivation that the presence of these off-diagonal terms derives from the fact that the 
background neutrinos are mass eigenstates, since they must be on shell. This origin of off-diagonal flavor space terms 
in the background potential due to neutrino- neutrino scattering was pointed out in Ref. |22|] . 

In calculations of the effects of neutrino flavor oscillations in the supernova environment, the supernova core is 
typically treated as a stationary source of neutrinos free-streaming from a "neutrinosphere," with the flux at the 
neutrinosphere being taken from large-scale numerical computations. Since the neutrinos forming the background 
also undergo flavor transformation, self-consistency between the oscillation probability and the background must 
be achieved. Following the framework of Ref. j23|, such a calculation was carried out in Ref. p3] in the quantum 
mechanical picture of neutrino oscillations. This involved a rather complicated procedure involving a flavor basis 
density matrix to describe the neutrinos above the neutrinosphere. Casual inspection of the form of the second term 
in Eq. (pTl) would seem to make this kind of approach necessary. However, having shown that this term can plausibly 
be written in terms of a macroscopic mass basis current, we see that the self-consistency between neutrino background 
and oscillation probability is most easily achieved by working in the mass eigenstate basis.^J Given effective interaction 
Hamiltonians like Eqs. ( p3|) a nd (^3|), the oscillation probability can be computed (in any basis) as described in Sees. 
FV| and 0. We write Eq. (g|) as 

f El dE p d( cos 9)d6 v ^ 
■W = J P j^y 3 [M(E P , cos 6, r) - f Pi (E p , cos 9, r)] |-, (46) 

where — (E p , pE p ), r is the radial position of a point above the neutrinosphere, and 9 is the angle between the 
neutrino momentum and the radial direction at the point with radial position r. Since the neutrinos are free-streaming, 
the distribution functions at r can be expressed simply in terms of the "known" neutrino distribution functions at 
the neutrinosphere, e.g.: 

f Vi {E p , cos 0, r) = f»j ( e p - cos 1>, R ) p - 3 -*« ( e p - cos 0> r ) . (47) 

where R is the radius of the neutrinosphere, and ip is the angle of the n eutrino emission w ith respect to the radial 
direction at the emission point; this angle is related to 9 by cosV- 1 = yl — [( r /P) sin#] 2 . The dependence of the 
oscillation probabilities on path length (and the background encountered on a particular path) are implicit in the 
r, 9 dependence. With an iterative procedure, self-consistency between the macroscopic neutrino currents and the 
oscillation probabilities should be achieved. 



It would seem reasonable to define neutrino flavor distribution functions in the relativistic limit. While at the emission point 
(i.e. the neutrinosphere) one could argue that these would be related to the mass basis distribution functions by f U( (p) = 
J"] \Uai | 2 fv a , at points above the neutrinosphere the relation between these sets of distribution functions is rather complicated, 
due to the flavor/mass oscillations of free-streaming neutrinos in a background. The resulting absence of a simple connection 
between the macroscopic flavor and mass neutrino currents at arbitrary position to plug into Eq. (^) makes working in the 
mass basis seem much more straightforward. 
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IV. CONSTANT BACKGROUND 



Having constructed effective interaction Hamiltonians as described in the last section, we employ the neutrino 
effective Lagrangian 



C = P[^(id fl -V IM P L )-M] u, 



(48) 



where M is the mass matrix and V M is a uniform background potential matrix. We make no assumptions about the 
number of neutrino generations or the structure of the potential matrix (other than to keep in mind that it might be 
singular). The canonical anticommutation relations yield the equation satisfied by the Green's function G(x,y), 



[7 M (^ - V„P L ) - M] G(x, y) = 5\x - y), 



(49) 



where i G(x, y) = (Tip(x)-ip(y))o. With our convention for the 7 matrices it is convenient to define the 2 x 2 (in spinor 
space) matrices G/j, the "chiral blocks" of the Green's function. Specifically Gjj is the nonzero 2x2 submatrix of 
PjGPj, where /, J can take the values L, R. 

In Sec. II we noted that, with the assumption of V-A interactions, Glr is the object of interest. We also saw in 
Sec. II that with the assumption of stationarity it is natural to Fourier transform the time variable while maintaining 
interest in the spatial dependence of the Green's function. Defining J = M~ 1 Grji, from Eq. ( |49"|) we find 



G L i{(w,x,y) = (u -Via ■ V)J(w,x,y), 



(50) 



where f(x,y) =J^e 



= f iHL p -iu(x°-y°) 



/(w,x, y). In the context of neutrino oscillation experiments we are in interested in 



well-separated source and detector positions, so we ignore terms in J with more than one factor of |x — y| in the 
denominator, 



d 3 p 
(27p 



(2tt) 3 , 

-iu|x— y I 



J(w,p) 



1 



w|x-y| 
J(uj,u,p = -f) - 



e H x - y| J(w,M,p = +f) 

J 



O 



u\x - y| 



(51) 



where we have integrated the cos# integral by parts, and defined u = |p| and f = (x — y)/|x — y|. It is evident that 
the two leading terms are azimuthally symmetric, and that their the sum is even in u. Furthermore, the Feynman 
boundary conditions should ensure that the two leading terms give equal contributions. We then have 



J(w,x,y) 



1 



/OO 
dwwe iu|x - y| J(w,w,p = +r). 
-00 



From Eq. 



(27r) 2 i|x-y| 

^), we find that J obeys the momentum space equation 

[uj 2 - |p| 2 - M 2 - loV° + p • V + cr • (V°p -LoV + iV x p)] J(lu, p) = 1, 



(52) 



(53) 



or D(uj,u,p = +f)J(w,u,p = +f) = 1. Detailed expressions for the spinor space elements of D(cu,u,p = +f) for 
general orientation of r are not particularly illuminating. However, it is easy to formally express the spinor space 
elements of J(w, u, p = +f ) in terms of the elements of D: 



J 1 V, u, p = +f) = [D 11 - D^iD 22 )-^ 21 ] ^ , 
J 22 (lo,u :P = +r) = [D 22 - D 2l {D ll )- l D 12 } ^ , 



J i2 (w,u,p = +f) = -(D 11 )-^ 12 ^ 2 



J (w,u,p 



22\-\jj21 jll 



(54) 
(55) 
(56) 
(57) 



These results are valid without any relativistic limit assumptions. Given specific mass and potential matrices, one 
could solve explicitly for J(u,u, p = +f). To study neutrino oscillations, we then need J(w,x, y), whose behavior is 
seen from Eq. (|5^) to be determined by the poles of J(w, u, p = +r ) with positive imaginary part (as determined by 
the Feynman boundary conditions). 
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A few general comments regarding these poles are in order. Consider for example J 22 , which can be expressed 

J 2 V, U ,p = +f) = (det D 11 ) [(det D ll )D 22 - D 21 [C ll ) T D 12 } _1 , (58) 



where (C ) is the transpose of the matrix of cofactors of D . Since the diagonal elements of D and D are 
second order in u, (det D 11 ) is of order 2n in u, where n is the number of neutrino generations; and overall the 
denominator of J 22 (cj,w, p = +f) will be a polynomial of order 4n in u. This is sensible in terms of a quasiparticle 
picture associated with the propagator: Each neutrino field, with two spin states each for particles and antiparticles, 
represents four states. For a single vacuum field the masses of these four states are degenerate; however, the presence 
of a parity and rotational invariance violating potential breaks this degeneracy. 

Let us examine the simplifications that occur in the relativistic limit and with source and detector localization. In 
the relativistic limit the poles contributing to the integral in Eq. (52) take the form 



2U 



(59) 



with the Feynman boundary conditions imposed by giving the "masses" rh 2 a small negative imaginary part. (There 
are also negative poles, with negative imaginary parts, that do not contribute to the integral; these factors each 
become ~ 2\u>\ when evaluated at the positive poles.) Furthermore, following the discussion of Sec. || regarding 
spatial localization, Eq. (|52|) takes the form 



•/(w,x,y) 



2|w| e^Mx-y) 
4tt|x s - y D | 



|xs— yo\ 




J [ u,u,p = rL 

M 



(60) 



where the sum is over the poles with positive imaginary parts. We recall that to is fixed by energy delta functions, 
to a positive value for neutrino oscillations and a negative value for antincutrino oscillations. Because the spatial 
localization sets f = ±L (with L taken to be the third spatial direction), the matrix D takes the relatively simple 
spinor space form 



D 



uJ,u,f>= — L 
M 



1 



r; ( yl 
— )v u - w 



iV 2 



1\V J 



(61) 



ur - u 2 - M 2 . 

Next we examine the momentum space pole structure of J 22 in the relativistic limit. In Eq. (|55|), since the residues 
of only the positive poles contribute, we can replace u by \uj\ whenever it multiplies a component of V 1 , committing 
errors of only O(uV /ui 2 ) or less with respect to other terms present. In that case, D 2 i vanishes and J 22 — » (D 22 ) -1 
with 



D 



22 



M 



-2q ■ V 



(62) 



where q = (|w|, L|w|) is the same as the neutrino momentum defined just below Eq. (^). We note that the denominator 
of J 22 is now only of order 2n in u; thus in the relativistic case, two of the quasiparticle propagating states are projected 
out. This is because in the relativistic limit the spin states naturally coincide with the chiral states. This is confirmed 
by noting that in Eq. (BBt), for example (w > case), as D 21 — * as u — > lu, half of the poles contributing to J 22 



come from (det D ) — > 0. But for 



D L 



— M 2 . Thus these poles correspond to the vacuum masses; 



these are the right handed particle states and left handed antiparticle states whose masses are unaffected by the left 
handed effective potential. As D 21 reaches zero, the contribution of these poles vanishes completely. 
For |w||x — y| 1, Eq. (|50"|) becomes 



G L i{(w,x,y) ~ w(l - a ■ L)J(w,x,y). 



(63) 



Since we have chosen L = (y^> — x<j)/|y£> — xg| to coincide with the third spatial dimension, the only nonzero 
component of Glr is G 22 R which for neutrino oscillations is G"j? R (\u>\;y,x) = 2\u>\ J 22 (|w|,y,x), where the spinor 
space indices are exhibited and the mass/flavor indices are suppressed (note that J 21 vanishes), and for antineutrino 
oscillations is G| 2 R (— |a>|; x, y) = — 2\uj\ J 22 (— \uj\, x, y). Thus we see that the Green's function takes the required form 
ofEqs. iflT 
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The matrix M 2 + 2q-V (or M 2 — 2q ■ V ) of Eq. (^2|) is precisely the effective mass matrix M 2 appearing in the usual 
quantum mechanical model of neutrino (or antineutrino) oscillations. The effective mass matrix can be diagonalized 
by a unitary transformation, UM 2 W = 1. Thus for neutrino oscillations, for example, 



(D 22 )^ = U^U* aj {u? -u 2 - rh 2 + ze)-\ (64) 
where Upj are the elements of U and rh 2 are the eigenvalues of M 2 . Then Eq. (|6^) becomes, using Eq. (|60|), 

e i|w|L-(y-x) 



G^(w,y,x) 



(65) 



where we have assumed that the various conditions discussed in Sec. II are satisfied. Comparison of Eq. ( |65| ) with 
Eqs. (17 24|) shows that the oscillation probability derived here is precisely the same as that found in the usual 
quantum mechanical model. The antineutrino case works out in a similar manner. 



V. NONUNIFORM BACKGROUND 



In this section we consider the case in which the effective potential — y M (x), that is, we allow it to vary in 
space (but not time). From Eqs. (|65|) , jl7|) , and (|24|), it is clear that in the constant potential case the portion of 
the Green's function comprising the oscillation amplitude obeys a Schrodinger-type equation, the same one used in 
the standard quantum mechanical picture. While one might think to simply replace the constant potential in this 
Schrodinger equation with a spatially varying one — thus arriving immediately at the standard result — we shall go 
back a little further in order to see what is being left out in the process. 

In allowing for spatial variation in V^, Eqs. (50-52J) are unchanged; but Eq. ( |53| ) becomes an integral equation, as 
J(w,p) must be convolved with the momentum space dependence of V 1 . Since such equations are difficult to deal 
with nonperturbatively, in this section, we take a different route of working with a partial differential equation in 
coordinate space. The coordinate space version of Eq. (l52|) is 



\lu 2 + V 2 - M 2 



;V°(x) -iV(x)- V-ier- (V°(x) V-iwV(x) + iV(x)x V )] J(w,x,y) = <5 3 (x-y) 



(66) 



Unlike in the case of an integral equation, to define the Green's function using this equation, we must also separately 
specify the boundary condition. We shall assume that the production region is localized to a region of adiabatically 
constant potential. Furthermore, since we expect the virtual particles to all have the same phase just after being 
produced, our boundary condition prescription will be that the the Green's function J asymptote to the constant 
potential Green's function on an infinitesimal sphere centered about yj^| 

The form of Eq. (B2J), together with our experience in the vacuum and constant potential cases, suggests that in 
pVf 2 /2|w| 2 <C 1, U|cl>|/|ci;| 2 <C 1 where we suppressed the matrix indices] we look for solutions of 



the relativistic limit 
J of the form 



With this substitution. 



l|x-y| 



47r|x — y 



■^(w.x.y). 



(V 2 



)J = <5 3 (x-y)e^ll*-ylF_ 
2U e iMI*-y' ■ 



VJ 



47r|x - 



2|o;| e *MI*-y| 
4?r|x - y| 
1 



-Lv 2 F + «(f.VF)-^ : 



(f • VF) 



iv 1 r 

—F-\ VF 

2 2w 2wx-y 



F 



(67) 

(68) 
(69) 



where as before f = (x — y)/|x — y|. Requiring the first term on the right hand side of Eq. (|6§| ) to cancel the delta 
function in Eq. (|66|) gives a boundary condition on F, namely (restoring flavor indices) 



5 0a . 



(70) 



3 Note that one must match more than just the limiting singularity of the Green's function at x = y to define a unique solution. 
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Aside from this boundary condition, we are interested in well-separated x and y (specifically |w||x— y| ^> 1), so that 
we may ignore the last term of Eqs. (J6^) and (|69|). Then Eq. ( p6| ) becomes 

.•(r.VF) + ^V^-^[M^ + ^- M (r.V) 

-*-(V°\w\i-uV + i\w\Vxi)]F + On^\VF\ \ = 0, (71) 

where in accordance with the relativistic condition I^|o;|/|a;| 2 <C 1, we will neglect the terms represented by 
C(V|a;| | VF|/|a;| 2 ) in comparison with the first term of Eq. (O). 

One can distinguish three cases: (1) |VF| ^> eF, where e = V° + M 2 /(2\u\) and M 2 denotes the largest mass 
matrix eigenvalue squared; (2) |VF| ~ eF; and (3) |VF| <C eF. One can argue that case (1) is not interesting since 
all terms leading to flavor mixing are rendered negligible. Case (3) is also not of present interest because one can 



argue using Eq. (71) that it violates our relativistic assumption. In case (2), |V 2 F/(2|oj|)| can be neglected compared 
with |VF|, provided that | W°\/e 2 < 1. Writing F' = r • VF, Eq. @ becomes 

iF' + ^-D(\u;\, X )F = Q, (72) 

where the spinor space elements of D(|cj|,x) for the particular case when x lies along the third spatial axis L arc 
given by Eq. (|6l| ) with V — > V(x) and u — ► \u)\. 

Since Eq. (p3f) holds under the assumptions of case (2), together with the spatial localization of the source and 



detector, we see that G 22 R is the only nonvanishing component. Hence, from Eqs. (|72|),(67),(p^|), and (B4J), we find 
that the neutrino oscillation amplitude H obeys the Schrodinger equation 

iH' = [M 2 + 2q ■ K(x)] H, (73) 
2|w| 

where q = (\lu\, L|w|). Similarly, in the case of antineutrino oscillations (lo < 0), the oscillation amplitude obeys 

iH' = ^ [M 2 - 2q ■ 7(x)] H. (74) 

Before concluding this section, a remark regarding the boundary conditions is in order. Note that F of Eq. fl72] ) 
satisfying the boundary condition Eq. (|7(]) is in general different from F satisf yin g Eq. (^36|) with its associated 
boundary condition [described just below Eq. (p6|)]. In particular, although Eq. ( |72| ) is valid naively only far away 
from y, as we threw out the last terms of Eqs. ( |68| ) and (^), we still insisted on the boundary condition Eq. ( fTOj ) at 
y to be the same as the boundary condition that would have been used for the exact equation. To justify this, we 
must show that the terms that we threw out are negligible even near the origin. We can argue this by noting that 
for case (2), we are already assuming |Vy°|/e 2 < 1 which turns out to imply (by expanding the potential to linear 
order in Taylor series about y) that in the relativistic limit the fractional variation of V° is much smaller than 1 until 
|w(x — y)| ^ 1 (after which the terms proportional to l/|x — y| that we threw out are negligible). That means that the 
potential can be treated as a constant until the terms proportional to l/|x — y| become negligible. This implies that 
one can place the boundary condition for the varying potential case on a sphere (centered about y) on which Eq. ( |7^ ) 
is valid using the solution to the constant potential case. As we saw in the last section, since the exact solution to the 
constant potential case on this sphere (with the a ppr opriate boundary condition) is, up to relativistically suppressed 
terms, the same as the solution obtained by Eq. (|7|) with Eq. (|70|), we can just set the boundary condition for the 
varying potential case using Eq. ( |70| ) as well. Note, however, that since our argument depends on the adiabaticity of 
the potential near the virtual particle production point, in other situations, one may need to be more cautious with 
the boundary conditions. 

Thus under appropriate conditions the results of the usual simplified picture are confirmed, including the boundary 
condition ^"(w, x, y)| y ^ x = 6?°. 

VI. CONCLUSION 

Starting from quantum field theory (QFT), we have defined a physically meaningful flavor oscillation probabil- 
ity; determined that portion of the neutrino propagator that comprises the oscillation amplitude; and derived the 
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"Schrodinger equation" for that amplitude in the presence of spatially varying background matter. As expected, 
the "Schrodinger equation" really corresponds to a time independent one since its derivation depends on the time- 
independence of the effective potential. In fact, the usual quantum mechanical approach is only really suited to 
problems in which oscillations occur in space only (that is, stationary systems like that studied here) or time only 
(e.g. the thermal bath in the early universe). While we have assumed the stationary case here, the basic framework 
could also be used to study oscillations in space in the presence of a time-dependent background. Ultimately, the 
description of flavor oscillating neutrinos in space and time in more general systems — i.e. those that do not lend 
themselves to interpretation in terms of a "source" and "detector" — would require a formulation in terms of density 
matrices (cf. Ref. p3|]) or Wigner functions (cf. Ref. [p^[). 

For situations that can be interpreted in terms of a "source" and "detector," we have also reviewed the conditions 
under which QFT will be useful in describing the neutrino oscillation process.]^] As long as we are looking in the regime 
in which the virtual neutrino goes on shell (rendering any propagator radiative corrections to be negligible or be merely 
a constant shift), the many body aspect of QFT is rendered irrelevant. In that case, only the production/detection 
vertex structure and the spins of the neutrinos are missing from the usual quantum mechanical treatment. These, in 
general, are less difficult to accommodate in the quantum mechanical treatment than the many body effects. Still, 
we believe them to be more straightforwardly accommodated in the quantum field theoretical treatment. For weak 
interactions, the chiral nature of the interactions combined with the relativistic nature of the on shell neutrinos 
suppresses all but one spin degree of freedom. Finally, the smallness of the neutrino mass splittings as well as the 
neutrinos going on shell allows one to factor out the production/detection part of the neutrino scattering process 
from the "oscillation" part, reducing the problem to the usual quantum mechanical system involving a single spinlcss 
particle. 

To state this another way, we have argued that in the context of stationary systems, the quantum field theoretic 
formulation used in this work is not of much use except when one or more of the following is true: the on-shell 
neutrino momentum is nonrelativistic, the production or detection vertices are non-chiral, the external particle wave 
packets vary appreciably about the pole value (value of the wave packet evaluated at the neutrino momenta) over 
momentum variations order of the inverse source-detector distance 1/L, or the effective mass splittings (determined 
by the effective potential including the background matter contributions) are large compared to the spread in the 
momenta of the external production/detection particles. 

When the neutrinos are non-relativistic or the interactions are not chiral, more than one spin contributes per 
amplitude. In that case the usual quantum mechanical treatment must be modified to incorporate the effects due to 
the various spin components. In particular, in this case an oscillation probability cannot be determined apart from 
the production and detection processes since two neutrino spin contributions are summed before the amplitude is 
squared. This is, of course, automatically accounted for in a quantum field theoretic treatment. Also, if the wave 
packet varies appreciably about the pole value when the momentum is varied by 0(1/L), the wave functions of the 
external particles creating and absorbing the virtual neutrino do not simply factorize out of the oscillation probability 
amplitude. Ref. |l5| is a study of one such situation. They find that the oscillation amplitude can exhibit a novel 
"plane-wave" behavior. Just as most effects, this probably can also be accounted for in the quantum mechanical 
formulation, but it is much easier in the quantum field theoretical treatment used in this paper. The involvement of 
the details of the external particles' wave function in determining the oscillation probability also applies when the 
effective mass splitting is much larger than the momentum spread in the external wave packets, but in such cases, 
it is not clear whether any neutrino oscillations can be observed in general because of the strong suppression of the 
amplitudes (however, see for example Ref. Jl5[|). 

We here make a few comments regarding Majorana vs. Dirac neutrinos. If the neutrinos were Majorana instead 
of Dirac, the only general arguments that would change are those dependent upon the existence of a right handed 
neutrino. Although we couched the mathematics in the Dirac spinor formalism, owing to the assumption of chiral 
nature of neutrino interactions and the relativistic limit, none of our general arguments depended upon the existence 
of a right handed neutrino. Hence, our conclusions for the relativistic limit are also valid for Majorana neutrinos. 
(Formulae for "neutrinos" apply to negative hclicity Majorana neutrinos, while those for "antineutrinos" apply to 
positive helicity Majorana neutrinos.) 

From our justification of the "Schrodinger equation" for stationary, flat spacetime systems, we expect that the 
heuristic ansatz used in Ref. pl| for studying neutrino oscillations in a stationary curved spacetime to be valid to 
the extent that the flat spacetime treatment is valid. In a nonstationary curved spacetime, in addition to the time 
dependence of the potential, there may arise extra complications of using the S-matrix formalism due to the nontrivial 



Many of these conditions were noted in Refs. pOHf 
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Bogoliubov transformations of the asymptotic states. This also deserves further investigation. 
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APPENDIX A: BOX QUANTIZATION AND LOCALIZATION 

In Sec. II we discussed the crucial role that wave packets play in the localization inherent in neutrino oscillation 
experiments. Studies of vacuum oscillations in quantum field theory typically begin with some kind of coordinate 
space description of the external particles, with the connection to the underlying momentum space wave packets 
not explicitly displayed. Since normalization in quantum field theory is generally fixed at the level of free-particle 
momentum eigenstates, fully normalized event rates cannot be computed unless the connection to momentum space 
is made. The usual practice in the literature has been to ignore this step, being satisfied with the identification of 
the factor called the "oscillation amplitude" in the simple quantum mechanical picture. The complete connection to 
momentum space is presumably a straightforward exercise, though perhaps tedious in the context of general wave 
packets. However, as shown in the text, the box wave packet allows a simple, reasonable approximation. For typical 
microscopic processes, a standard trick used to relate the S matrix to rates and cross sections is to employ the 
fiction of a box with quantized states, which reproduces the box wave packet results. For pedagogical purposes, we 
here demonstrate that the use of a similar fiction — separate boxes confining the external particles at the source and 
detector — leads cleanly and directly to an event rate in the form of Eq. (|l]) . In this heuristic approach the boxes serve 
the double duty of imposing the localization of source and detector necessary to the existence of oscillations, and the 
normalization that gives an absolute event rate. 

For illustration we will consider a particular neutrino oscillation process. The first vertex of the process we consider 
involves the scattering of a charged lepton of flavor a off a proton in a source region of volume Vs, producing a final 
state neutron confined to the source region and a virtual neutrino that propagates over all space. At the second 
vertex, the virtual neutrino interacts with a neutron confined to a detection region of volume Vb, producing a proton 
and charged lepton of flavor (3 confined to the detection region. The free Lagrangian is taken to be 



V / d 3 x^ s (i-f t "d^ l - m)4> s + V / (fxi/jdii-y^dfj, - m)ip d 

s JVs d JV D 

+ d 3 xMird^-m)i> t , (Al) 



Lo = 



all space 

where the index s runs over the particles confined to the source (a,p, n), the index d runs over the particles confined 
to the detector (/3,p, n), and the index i runs over the neutrino mass eigenstates. The box localization is really a result 
of complicated classical potentials preparing the initial state. The expansions of the source and detector particle fields 
ard3 

1>.A X ) = -/fr= -/fF= \ a P,*,(s,d)U pMs ,d)e- ip - x + ^, (M) ?W,(M) elP i > ( A2 ) 



/2E P 

v - ' - (j p V ' 

where the momentum sum is over the discrete momenta arising from the application of periodic boundary conditions 
to the box, the index a is over the spin states, and the commutation relations of the creation/annihilation operators 
are of the form [a Py(T , s , a p , a , s ,] = <5 PP ' 6 a <j' 5 S . S >. In contrast, the free neutrino fields have a continuous momentum 
spectrum, 



Since we are taking momentum quantization "seriously" in this approach, it is convenient to adopt slightly different nor- 
malization conventions from the ones used in the main text. The differences are easily deduced from the explicit expressions 
for the field expansions and the commutation relations given above. 
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/d 3 v 1 
(2tt) 3 / 2 CT ' *^ P ' CT ' *) e ~*' S + 6t ( p ' CT ' *M p ' CT ' O^""] ' 

with [a(-p,<T,i),at(p',a',i')] = S 3 (p - p') 8 aa , 8 iV . 

The interaction Lagrangian relevant to the first vertex of the process described above is 



(A3) 



Lj = -A / d X%l) a ^(l ~ 7 5 )^aV'p7/i( 1 - 9Al5)lpn + H.< 



(A4) 



where g = Gf cos 9 c /y/2, Gf is the Fermi constant, 9 C is the Cabibbo angle, and g A ~ 1-26 is the neutron-proton 
axial vector coupling. A similar interaction Lagrangian describes the interaction in the detector. The relation between 
the flavor fields and mass eigenstate fields is v a = U a iipi, where the U a i are elements of a unitary matrix. The 
amplitude for the neutrino production/detection process is 



T a p = -g 



-ip-x 



dx° / dy° / d 3 x / d 3 y- ^_ _ 



-il-y 



oil' -V 



cik'-y 



[fi(lV(l - 9AJ5)U(1)] 



VVhVm VVEVzW VvEV^ 

x [s(k')7 Al (l-75)^ Q (y,x)) 7l ,(l-75)«(k)] [u(p')^(l-g Al5 )u(p)}. 



(A5) 



The propagator or Green's function, iG^ a (y,x) = (T{vp(y)i/ a (x)})o, with T{} and (}o denoting a time-ordered 
product and vacuum expectation value respectively, can be expressed in momentum space as 



d 4 s 



-is-(y-x) Q0<xf s \ 



(A6) 



Integration over a: , y°, and s° in Eq. ( |A5| ) fixes \s°\ = \E k + E p — E p >\ and yields a factor (2ir) 2 5(E\ i + E p + E\ — 
Ek' — Ep/ — E\> ) . 

We let Vs and Vd be cubes of sides L$ and Lp, centered on and yo respectively, and employ the approximation 
of Eq. (fjl|). Because of the V-A lepton currents, the relevant block of the Green's function is Glr- Suppose that 
this block of the propagator takes the form of Eq. ( |l7|), with H^ a having only flavor indices. Changing variables to 
x' = x — xg and y' = y — yo , the amplitude in Eq. (|A5|) can be expressed 



T, 



a/3 



i\ 2 6(E k + E P + E Y - Ey - E p , - Ey) 

2{v s f/ 2 {v d f/ 2 \y D - Xs | (n s Vie;) (u d vm) 



e i(k+p-p')-xs e i(l-l'-k')-yD 



:^(£ q Jfl,x S ) 



Y[A(u I ,L s )A(v I ,L D ) 



[u(p')^(l-g Al5 )u(p)] 



x [u(k') 7ti (l- l5 )u(q)-] [^(q) 7 ,(l-7 5 )u(k)] [u(l') 7 »(l - g Al5 )u(l)} 



(A7) 



where the index I is over the three momentum directions, u = k + p — p' — q, v = 1 — 1' — k'+q, A(w, a) = 
(2/w) sin(uia/2), and q = Eq L. 

An event rate is obtained by squaring the amplitude; interpreting one of the energy delta functions as T/(2n), 
where T is a time interval; and dividing by T . At this stage we will also make the approximation of a continuum 
of states in the source and detector, i.e. take (Lg/2), (Lu/2) to be large. In this limit A(w,a) — > (2ir)5(w), and 
[A(w, a)] 2 — > [(2/w) sin(wa/2)| to ^o] (2ir)5(w) — 2na S(w). Finally, in the continuum limit there is a phase space factor 
of the form d 3 pV/(2Tr) 3 for each final state particle. The event rate obtained from Eq. (A7) is then 



rfr 



a/3 



A 4 6(Ek + E P + E Y - Ev - E p , - Ey,) 



A{2^V S L 2 (Jl s 2E s ){Y{ d 2E d ) 
x<5 3 (p + k p' q)<5 3 (l + q - 1' - k') d 3 p' d 3 l' d 3 k' 

x [|S(p')7 M (l-ffA75)^(p)S(k')7^(l-75)«(q)| 2 



spins 



x |u(q) 7 ,(l - 75Mk)u(l')7' y (l - ffA75)u(l)r 



(A8) 
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Using standard plane wave methods to compute the neutrino production rate and cross section factors in Eq. ([I]), it 
is easy to verify that Eq. (AS) is equivalent to Eq. ([!]), with P v<x ^ v& — \H^ a (E <i , yjj, xs)| . As noted in the text, 



explicit computations yield the oscillation probabilities found in the usual quantum mechanical model. 



Eq. ( |A8| ) represents the rate per source proton, source charged lepton, and detector proton. To get the total 
experimental rate it is necessary to sum over these source and detector particles. Employing the usual classical 
distribution function, for uniform spatial distribution the number of particles of type x is 

N * =v Jwf f * {p) > (A9) 

so that the total event rate is 

r rfi v (pk d 3 l 

dT cxp = J V S j^ f P (p) V S j^ / Q (k) Vn-^ f P (l) dT aP . (A10) 

We note that one of the factors of Vs cancels, leaving one factor each of Vs and Vd- These remaining two volume 
factors are replaced by d 3 xg and (Pyo, and integrations over these source and detector positions are performed. While 
we have explicitly shown how this works out for our particular chosen process, for any neutrino production/detection 
process one factor of Vs and Vd will remain at the end. 
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